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1.0 


INTRODUCTION 


1.1.  TASK  I 

The  ionic  liquid  (IL)  program  at  the  Air  Force  Research  Laboratory  is  investigating  tailored  energy- 
dense  liquids  to  provide  a  scientific  foundation  for  the  advancement  of  the  performance  and 
operability  envelopes  of  current  propulsion  systems.  General  efforts  involve  the  discovery  of 
energetic  ILs  based  on  heterocyclic  and  open-chain  cations  in  combination  with  reactive  anions.  A 
major  goal  has  been  the  development  of  IL  fuels  which  undergo  hypergolic  ignition  upon  contact  with 
common  propulsion  oxidizers.  In  state-of-the-art  bipropellant  systems,  N2H4  is  the  fuel  and  N204  is 
the  oxidizer.  Both  compounds  are  highly  toxic,  hydrazine  being  a  suspected  carcinogen.  Replacing 
hydrazine  with  an  ionic  liquid  which  generally  possesses  very  low  volatility,  and  therefore  virtually  no 
vapor  toxicity,  would  be  extremely  desirable  if  not  already  critical.  At  the  same  time,  the  toxic  oxidizer 
should  be  exchanged  for  a  high-performing,  environmentally  benign  oxidizer  like,  perhaps,  hydrogen 
peroxide  (decomposition  of  hydrogen  peroxide  affords  only  oxygen  and  water),  if  other  desired 
properties,  especially  hypergolic  reactivity  with  the  fuel,  can  be  retained.  So  far,  most  researchers 
seeking  hypergolic  fuels  have  limited  themselves  to  the  extremely  toxic  and  corrosive  nitric  acid 
solutions.  While  important  questions  remain  unanswered,  we  are  exploring  new  ground.  During  our 
previous  work  we  demonstrated  the  hypergolicity  of  ILs  with  one  class  of  hydrogen-rich  anions 
toward  a  variety  of  common  propellant  oxidizers,  including  hydrogen  peroxide.  However,  for 
practical  purposes  the  materials  had  several  shortcomings,  especially  poor  liquid  range  and  viscosity. 
Furthermore,  their  high  hydrocarbon  content  limits  propellant  performance.  Based  on  multiple 
molecular  design  strategies  the  proposed  research  effort  focuses  on  synthesizing  novel  energetic  ILs, 
supported  by  experimental  mechanistic  studies  and  computational  investigations.  Expanding  on  our 
previous  work  we  intend  to  enlarge  the  chemical  landscape  of  known  light  metal  hydride  ILs.  We 
envision  the  successful  completion  of  this  project  will  result  in: 

■  Fundamental  understanding  of  the  coordination  of  light  metals  in  highly  reducing 
environments. 

■  Expansion  of  the  range  of  materials  available  for  hydrogen  storage. 

■  Novel  solvent  systems  enabling  unique  reduction  reactions. 


1.2  TASK  II 

The  Ionic  Liquids  Ignition  program  investigates  the  chemical  kinetics  and  reaction  dynamics 
involved  in  the  hypergolic  and  catalytic  ignition  of  ionic  liquid  propellants  with  the  purpose  of 
identifying  key  intermediates  and  kinetic  bottlenecks  which  can  enhance  or  restrict  performance  in 
ionic  liquid-based  propulsion  systems.  In  better  understanding  the  energy  landscape  involved  in 
hypergolic  ignition  of  ionic  liquids  with  various  oxidizers,  we  can  use  this  knowledge  to  support  the 
synthetic  efforts  in  Task  I  to  discover  viable  candidate  ionic  liquid  fuels  for  propulsion  systems.  An 
important  goal  of  our  research  is  to  use  sensitive  and  selective  experimental  probes  to  understand  in 
real  time  and  at  the  molecular  level  the  underlying  chemistry  involved  under  relevant  extreme 
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environments.  The  experimental  data  will  provide  the  necessary  mechanistic  and  chemical  kinetics 
information,  and  with  ab  initio  quantum  chemical  analysis,  accurate  reaction  mechanisms  will  be 
elucidated.  The  combined  information  will  feed  into  chemical  kinetics  models,  which  will  make 
reactivity  predictions  possible  for  hypergolic  fuels  and  high-temperature  catalysts,  leading  to 
improved  design  of  energetic  ILs  including  more  reliable  ignition  and  sustained  combustion. 


2.0  TECHNICAL  SUMMARY 
2.1  TASK  I 

Propulsion  performance  can  be  fostered  by  light  metals  which  have  large  combustion  energies 
and  relatively  light  products.  Elements  with  considerable  performance  advantages  and  non-toxic 
products  are  aluminum  and  boron.  In  addition,  high  hydrogen  content  fulfills  the  need  for  light 
combustion  products  through  the  production  of  hydrogen  gas  and  water  vapor.  ILs  containing 
AI(BH4)4  anions  may  be  viewed  as  a  densified  form  of  hydrogen  stabilized  by  metal  atoms.  The 
volumetric  hydrogen  content  of  tetraethyl-  ammonium-AI(BH4)4  is,  in  fact,  99%  higher  than  that 
of  liquid  hydrogen.1  Previously  we  reported  an  IL  containing  the  [AI(BH4)4]'  anion  but  the  material 
has  several  shortcomings  which  make  it  impractical  to  use  as  a  propellant.  In  an  attempt  to  lower 
the  viscosity,  neutral  ABH  was  mixed  with  an  IL  containing  the  cyanoborohydride  anion.  ILs  with  the 
[NCBH3]'  anion  by  itself  generally  possess  very  low  viscosities  (~20cp).  It  was  anticipated  that  the 
combination  of  [NCBH3]"  with  ABH  would  produce  a  new  anion  of  the  formula  [AI(BH4)3NCBHb]\  A 
couple  of  possibilities  were  considered  for  the  binding/coordination  of  the  [NCBHs]"  anion  to  ABH. 
Boron  and  aluminum  can  compete  for  the  carbon  and  nitrogen  and  calculations  revealed  that  boron 
prefers  carbon  by  ~13kcal/mol.  Direct  coordination  of  the  CN  group  to  aluminum  is  also  preferred 
over  3-centered-2-electron  hydrogen  bonds  by  ~21kcal/mol  (Figure  1). 

However,  during  our  previous  studies  we  discovered  that  the  reaction  chemistry  of  ABH  is  far  more 
complex  than  initially  anticipated  and  our  initial  findings  prompted  us  to  investigate  this  chemistry  in 
more  detail. 


Figure  1.  Different  coordination  possibilities  of  [NCBH 3]'  to  ABH. 

We  established  the  existence  of  a  prevalent  mixed  borohydride/cyanoborohydride  aluminum 
compound  with  a  coordination  number  of  seven  (Figure  2)  and  saw  spectroscopic  evidence  that  other 
species  are  present. 
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Figure  2.  A  new,  doubly  charged  anion  (calc.,  left;  crystal  structure,  right)  [AI(BH4)2(NCBH3)3]2' 

In  a  next  step  we  first  explored  the  synthetic  accessibility  of  new  anions  of  aluminum  by  reacting 
ABH  with  a  large  amount  of  [NCBH3]'. 


Scheme  1.  Reaction  of  ABH  with  Methyl-triphenyl-phosphonium  CBH. 

Final  composition  of  the  only  isolated  product  was  resolved  by  a  single  crystal  X-ray  structure 
determination  and  the  structure  of  the  novel  [AI(CNBH3)6]3‘  anion  is  shown  in  Figure  3. 


Figure  3.  X-ray  crystal  structure  of  isolated  product  showing  a  rare  example  of  a  hexacoordinated, 
triply  charged  aluminum  anion. 
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In  2012  the  group  of  Ingo  Krossing  investigated  literature  known  procedures  to  make  heterocyclic 
borohydride  salts2  in  great  detail  and  determined  that  the  best  material  obtainable  thereby  still  had 
a  halide  content  of  at  least  22. 5%.3  In  addition  they  discovered  a  new  method  for  preparing 
analytically  pure  borohydride  salts  in  a  mixed  solvent  system,  liquid  ammonia/methylene  chloride,  at 
low  temperatures. 

We  now  discovered  another  method  which  allows  us  to  run  the  metathesis  in  a  single  solvent  at 
room  temperature. 

We  took  advantage  of  the  unique  coordination  possibilities  of  alkali  metal  borohydrides.  The 
coordination  complex  between  methylimidazole  and  sodium  borohydride  possesses  a  distinct 
characteristic  which  is  not  found  in  uncoordinated  NaBH4  (Scheme  2). 


r=\ 


H3C— n\^n 


+NaBH4 


r.  t.  4  days 


-> 


f=\ 

-N^N 


NaBH4 


Scheme  2.  Coordination  chemistry  of  methylimidazole  with  NaBH4. 


While  NaBH4  is  soluble  only  up  to  ~2%  in  acetonitrile  the  coordination  complex  with 
methylimidazole  revealed  a  high  solubility  in  this  solvent.  Many  heterocyclic  halide  salts  are  soluble 
in  acetonitrile  as  well.  Therefore,  the  new  coordination  complex  could  be  used  in  metathesis  reaction 
with  the  formation  of  sodium  chloride  and  the  free  amine  (Scheme  3).  The  free  amine  could  be  easily 
removed  by  washing  the  product  with  diethylether. 


f=\ 

-M^M 


r=\ 

NaBH4  +  H3C — N  — C2H5 

>  I  cr 


ch3cn 


12-16  h,  r.t. 


r\ 

H3C — N  .N— C2H5 

© 


BH/ 


CH3  ch3 

Scheme  3.  Metathesis  reaction  between  heterocyclic  halide  salt  and  NaBH4  coordination  complex. 


Our  new  method  has  multiple  advantages  over  the  procedure  developed  by  Krossing.  Krossing  has 
to  rely  on  a  co-solvent  system,  between  ammonia  and  methylene  chloride.  Reactions  have  to  be 
carried  out  at  low  temperature.  Methylene  chloride  has  more  handling  restrictions  and  health  issues 
than  acetonitrile.  This  new  synthetic  route  to  analytically  pure  heterocyclic  borohydride  salts 
provides  the  necessary  starting  materials  for  an  easy  conversion  to  a  new  class  of  ILs  with  complexed 
[AI(BH4)4]‘  anion. 

Though  the  melting  points  (MP's)  of  the  borohydride  derivatives  were  high  (Scheme  4)  it  was 
anticipated  that  the  liquefying  nature  of  the  aluminum  borohydride  anion  should  allow  for  a 
significant  MP  depression.  Upon  complexation  two  of  the  three  salts  formed  viscous  room 
temperature  ionic  liquids  (RTIL's)  (Scheme  4,  compound  1  and  2),  while  3  with  only  a  44°C  MP 
depression  was  not  a  RTIL. 
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l-ethyl-2,3-dimethylimidazolium  tetrakis-tetrahydroborate  aluminate  (3)  was  dissolved  in  a  minimal 
amount  of  dimethoxyethane  (DME)  to  carry  out  drop-test  experiments.  While  the  initial  tests  were 
successful  a  precipitate  was  formed  after  a  couple  days.  The  precipitate  consisted  of  a  new  complex 
in  which  two  borohydride  ligands  were  displaced  and  instead  a  bidentate-p2-bridging 
methoxyethoxide  ligand  system  was  installed.  This  is  the  first  known  example  of  such  a  mixed  ligand 
system  (Scheme  5). 


In  order  to  establish  that  this  unusual  coordination  is  general  for  aluminum  borohydride  containing 
materials,  the  etherate  of  aluminum  borohydride  was  reacted  with  DME  (Scheme  5).  The  material 
recovered  from  this  reaction  proved  to  be  identical  to  that  from  the  previous  salt  solution.  This 
demonstrates  that  aluminum  borohydride  containing  compounds  are  not  compatible  with  strongly 
coordinating  solvents  like  DME.  This  will  have  to  be  considered  for  future  development  of 
salt/solvent  systems  which  could  deliver  useful,  high-hydrogen  materials. 
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Scheme  5.  Reaction  between  Aluminum  borohydride-etherate  and  DME. 
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Computational 


Preliminary  simulations  have  uncovered  some  of  the  rich  chemistry  inherent  in  the  aluminum 
borohydrides.  Oxidation  by  water  of  either  the  simple  borohydride  anion  or  the  neutral  borane 
molecule  in  the  gas  phase  produces  hydrogen  gas  with  the  formation  of  a  boron-oxygen  bond  and 
takes  place  via  the  transition  states  displayed  in  Figure  4. 


Figure  4.  Transition  states  of  BH4-H2O  and  BH3-H2O. 


It  appears  that  oxidation  of  the  aluminum  borohyride  anion,  [AI(BH4)4]',  does  not  occur  directly  in  an 
analogous  manner.  Rather,  a  two-step  reaction  in  which  water  initially  displaces  a  neutral  borane 
molecule  (Figure  5,  left)  and  subsequently  reacts  with  it  in  a  complex  with  the  remaining  anion  (Figure 
5,  right)  is  observed  in  the  calculation. 


Figure  5.  Energy  profiles  along  the  Intrinsic  Reaction  Coordinate  (IRC)  for  the  hydrolysis  of  the  [AI(BH4)4f. 

Additional  insight  into  the  relative  reactivities  of  these  species  is  gained  from  the  calculated  activation 
barriers  and  reaction  energies  given  in  Table  1.  Although  hydrolysis  of  [AI(BH4)4]'  nominally  involves 
oxidation  of  BH3,  the  presence  of  the  intermediate  anionic  fragment  leads  to  an  activation  barrier 
and  reaction  energy  in  between  those  for  bare  neutral  BH3  and  the  bare  [BH4]'  anion.  These  values 
agree  with  the  facts  that  neutral  boranes  are  pyrophoric  in  moist  air,  that  many  aqueous  borohydride 
solutions  have  substantial  stability,  and  that  aluminum  borohydride  ILs  manifest  considerable 
hydrolytic  stability,  perhaps  intermediate  between  the  two.  Theory  will  help  the  experimental  effort 
by  predicting  early  intermediates  such  as  [H-AI-(BH4)3]',  and  modeling  their  vibrational  spectra  as  an 
aid  to  identifying  the  species  evolving  in  cryogenic  matrices. 
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Table  1.  Energies  (kcal/mol)  for  the  reactions  which  form  H2  by  hydrolysis. 


Activation  energy  Ea 

Reaction  energy  AE 

BH3 

39 

-21 

[H-AI-(BH4)b]~  +  bh3 

46 

-12 

[BH4]“ 

69 

+4 

TASK  I  Specific  research  goals  and  objectives  for  FY16. 

Unfortunately,  neither  the  cyanoborohydride  approach  nor  the  pure  tetrakis-tetrahydroborate 
aluminate  ILs  discovered  possess  the  desirable  physical  properties  for  a  suitable  propellant.  In 
general,  their  major  shortcomings  still  are  poor  liquid  range  and  high  viscosity  as  well  as  the  high 
hydrocarbon  content  of  the  cations  required  to  lower  their  melting  points  which  severely  limits 
propellant  performance. 

In  light  of  our  discoveries  we  turned  to  systems  simpler  and  higher  performing  than  the  complexed 
AI(BH4)4‘  anions.  Solutions  of  Li-AI  hydrides  and  LiBhU  in  ethers  have  shown  some  desirable  propellant 
properties,4  and  a  variety  of  lithium  metal  hydrides  are  commercially  available.  This  suggested  an 
entry  point  into  the  coordination  chemistry  of  lithium  salts  with  some  new  high  energy  heterocyclic 
ring  systems.  Our  studies  on  some  of  these  promising  new  heterocyclic  systems  is  summarized  below, 
however,  the  objective  for  this  last  performance  period  were: 

1.  Study  stable  lithium  borohydride  tetrazole  complexes. 

2.  Study  lithium  borohydride  complexes  of  simple  amines  with  high  borohydride-to-amine  ratios. 


Synthesis  of  5-substituted  tetrazoles  as  novel  ligands  and  energetic  salts 

In  pursuit  of  functionalized  energetic  heterocycles  suitable  to  form  stable  coordination  compounds 
with  metal  borohydrides,  the  highly  energetic  tetrazole  ring  system  was  chosen  and  it  was  desired  to 
incorporate  different  hydrazine  moieties.  Hydrazines  are  still  the  State  of  the  Art  propellants  and  it 
was  hoped  to  retain  some  of  their  desirable  properties  in  a  new  hydrazinotetrazole  derivative. 


At  first  it  was  envisioned  that  simple  nucleophilic  substitutions  would  be  a  straight  forward  approach 
to  prepare  5-(hydrazino-alkly)  tetrazoles  (Scheme  6). 
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Scheme  6.  Synthesis  of  5-(hydrazino-propyl)  1H  tetrazole  3  and  trimethylene  tetrazole  6 . 
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Disappointedly  our  initial  approach  of  a  single  step  replacement  of  the  halide  group  with  a  hydrazine 
moiety  leading  to  pure  5-(hydrazino-alky)-lH-tetrazoles  1  and  2  were  not  successful.  Thus  an 
alternative  multistep  synthetic  route  involving  the  use  of  BOC  protected  hydrazine  as  a  substrate  was 
devised  (Scheme  7). 


Scheme  7.  Synthesis  of  hydrochloride  salts  9  and  10 

Initial  attempts  to  neutralize  the  hydrochloride  salt  9  with  standard  alkaline  bases  as  well  as  metal 
alkoxides  led  to  inseparable  mixtures  of  inorganic  salt  and  the  desired  5-(hydrazino-methyl)  tetrazole 
1.  A  lengthy  washing  and  filtering  procedure  was  employed  which  finally  afforded  5-(fiydrazino- 
methyl)-l-H-tetrazoles  1  as  clean  white  powder  (Scheme  8).  X-ray  quality  crystals  of  1  were  obtained 
by  diffusion  of  diethyl  ether  into  methanol  solution  (Figure  6). 
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Scheme  8.  Synthesis  of  S-fhydrazino-methyl),  5-  (hydrazino-ethyl)  tetrazoles  1  and  2 


Figure  6.  Single  crystal  X-ray  structure  of  the  hydrochloride  9  and  5-(hydrazino-all<ly)  tetrazoles  1. 
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The  coordination  chemistry  of  lithium  borohydride  has  been  studied  previously.  Almost  70  years 
ago  its  coordination  with  ammonia5  was  investigated  and  more  recently  much  work  was  carried  out 
emphasizing  the  structural  chemistry  of  lithium  borohydride  with  various  types  of  coordinating 
ligands.6'9  Most  notably  Noth  et.  al.  have  developed  a  methodology  for  describing  the  interaction  of 
the  hydrogens  with  the  metal  center  in  coordination  compounds.  His  work  demonstrated  the  very 
complicated  relationship  between  the  borohydride  and  the  metal,  influenced  by  the  electronic  and 
steric  demands  of  the  ligands  (Figure  7). 6 


Figure  7.  Modes  of  borohydride  coordination  with  metal  centers  in  solvated  systems 


Our  preliminary  ventures  into  the  coordination  chemistry  of  lithium  borohydride  have 
concentrated  on  the  coordination  of  high  nitrogen  ligands.  Over  the  last  decade  high  nitrogen 
compounds  have  become  fundamental  to  the  search  for  new  energetic  materials. 

1,5-dimethyltetrazole  forms  a  complex  with  a  1:2  lithium  to  ligand  ratio  (Scheme  9).  It  is  interesting 
to  note  that  the  ligand  bridge  the  two  lithium  metal  centers  (Figure  8).  This  is  the  first  example  where 
two  different  nitrogen  atoms  of  the  same  ring  coordinate  and  bridge  multiple  lithium  cations. 
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Figure  8.  Crystal  Structure  of  Li  (1.5-dimethyltetrazole)2  BH4 


Interestingly,  a  modest  change  in  the  steric  interactions  by  isomerization  to  the  2,5- 
dimethyltetrazole  derivative  resulted  in  a  1:1  complex  (lithium  to  ligand)  (Scheme  10). 
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Scheme  10.  Coordination  of  2,5-dimethyitetrazoie  with  lithium  borohydride 
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The  resulting  polymeric  structure  shows  extended  lithium  and  borohydride  interactions.  To  our 
knowledge  this  is  the  first  example  of  this  type  of  bonding  motif  with  lithium  borohydride  (Figure  9). 


Figure  9.  Crystal  structure  of  Li  (2,5-Dimethyltetrazole)BH4 

In  the  context  of  the  work  by  Noth  et.  al.,  the  two  different  coordination  complexes  confirm  the 
very  diverse  role  that  the  borohydride  ligand  plays.  More  important  from  a  propellants  perspective 
is  the  change  in  stoichiometry  of  the  performance  enhancing  borohydride  with  such  a  small  change 
in  structure.  Although  both  complexes  are  solids  at  ambient  temperature,  our  overall  goal  is  to 
solvate  LiBFU  at  a  high  mole  ratio  and  the  general  presence  of  multiple  bonding  modes  might  aid  in 
liquefying  these  systems.  As  a  next  step  coordination  with  the  low  boiling  amine,  methylamine,  was 
investigated.  A  2:3  lithium  to  amine  ratio  was  determined  by  1H  NMR  (Figure  10).  No  undesirable 
amino  borane  side  product  was  observed  as  can  be  seen  by  the  single  pentet  present  in  the  nB  NMR 
(Figure  11). 


-30  -40  -50 


90  80  70 


Figure  11.  11 B  NMR  of  Li2(H2NCH3)3(BH4)2 
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Surprisingly,  this  complex  showed  a  melting  point  of  only  4°C  and  appeared  to  be  thermally  stable 
with  no  decomposition  visible  until  ~180°C  according  to  DSC  data  (Figure  12).  Upon  further  analysis 
by  TGA-MS,  it  was  seen  that  the  coordination  of  the  amine  was  not  as  strong  as  initially  supposed. 
The  methylamine  ligand  (mass  31)  is  almost  immediately  lost  upon  heating  (Figure  13).  Flowever, 
this  represents  the  first  example  of  a  room  temperature  liquid  LiBFU  complex. 


Figure  12.  DSC  of  Li2(H2NCH3)3(BH4)2  Figure  13.  TGA-MS  of  U2(H2NCH3)3(BH4)2 


At  this  stage  in  our  preliminary  studies  we  decided  to  investigate  other  ligands  with  chelating 
abilities.  The  first  ligand  investigated  was  diethylenetriamine  as  an  analogue  to  the  known  Li-glyme 
complexes  (Figure  14). 


Figure  14.  Solvation  by  diethylenetriamine  compared  to  the  known  Li-glyme-complex. 


According  to  NMR  spectroscopy  a  1:1  lithium  to  ligand  complex  was  formed  (Figure  15). 
Observations  with  water  as  the  solvent  for  the  initial  NMR  were  tantalizing.  While  NaBFU  solutions 
are  known  to  be  water  stable  for  quite  some  time  at  a  pH  of  10,  LiBFU  has  marginal  hydrolytic  stability. 
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Even  though  the  water  used  for  our  NMR  samples  had  a  pH  of  4,  we  were  surprised  to  note  very  little 
gas  evolution  upon  dissolving  the  complex  and  the  immediately  recorded  1H  and  nB  NMR  spectra 
showed  only  very  slight  hydrolysis.  Further  spectra  taken  after  2  days  and  5  days  showed  only 
minimal  hydrolysis  of  the  UBH4  complex  (Figure  16).  We  believe  this  is  the  first  time  a  water  stable 
complex  of  LiBH4  has  been  observed  (at  least  short  term)  and  that  these  complexes  should  certainly 
be  investigated  for  their  ability  to  be  practical  reducing  agents  in  water  solutions. 


Figure  15.  Initial  1H  (left)  and 11 B  NMR  (right)  of  Li(diethylenetriamine)BH4 


i _ j 

_ 1 _  1  1 

- 

16  15  14  13  12  11 

10  987654321  0  -1  ppm  90  80  70  60  50  40  30  20  10  0  -10  -20  -30  -40  -50 

y  Sr  iii  w 

ll  1 

i . JL . _ . . . ini .  ,  i 

-60  -70  -80  ppm 

5.0  4.5  4.0  3.5  3.0  2.5  2.0  1.5  1.0  0.5  0.0  -0.5  -1.0  -1.5  ppm  90  80  70  60  50  40  30  20  10  0  '10  "20  '30  '40  '50  '60  "70  '80  PPm 

"'iirJ  III  V 

Figure  16. 1H  (left)  and 11 B  NMR  (right)  of  Li(diethylenetriamine)BH4  after  2d  (bottom),  after  5d  (top) 
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Collaborations 


An  important  collaboration  was  established  with  Purdue  University  (Prof.  Timothee  Pourpoint). 
Under  an  AFRL-funded  SBIR  contract  a  small  thrust  stand  was  developed  by  Orbitec  (Orbital 
Technologies  Corporation)  and  delivered  to  Purdue  University.  This  government-furnished 
equipment  (Figure  17)  was  modified  during  a  2015  effort  funded  by  our  laboratory  to  allow  it  to 
accommodate  our  new  fuels  and  complementary  oxidizers.  This  test  stand  will  be  used  to  learn  more 
about  the  combustion  characteristics  of  molecular  fuels  containing  light  metals  and  will  therefore 
quantify  how  much  of  the  promised  performance  is  likely  to  be  possible.  This  constitutes  an  essential 
step  in  guiding  the  transition  from  research  synthesis  to  propellant  development.  The  miniature  test¬ 
bed  thruster  will  be  operated  at  about  1  to  5  Ibf  with  an  impinging  jet  configuration.  The  purpose  of 
this  thruster  is  to  allow  for  performance  characterization  in  terms  of  C*  and  Isp.  It  has  a  flexible 
design  that  permits  the  use  of  replaceable  injectors,  nozzles,  and  combustion  chambers  of  various 
sizes.  Measurements  include  thrust,  chamber  pressure,  and  chamber  temperature.  Some  of  the 
preliminary  2015  test  results  are  depicted  below  showing  a  drop  test  result  (Figure  18)  and  successful 
test  firing  of  a  boron  containing  fuel  developed  during  our  efforts  (Figure  19  left).  The  research  effort 
will  include  optimizing  system  parameters  as  well  as  oxidizer  to  fuel  ratios  to  demonstrate  acceptable 
combustion  efficiencies  not  realized  during  our  initial  tests  (Figure  19  right). 


Figure  17.  Modified  test  stand  installed  at  Purdue  University,  Zucrow  Test  Facility 


Figure  18.  Drop  test  image  ofAFRLfuel  and  WFNA:  Ignition  withion  1.1ms 
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Figure  19.  Successful  hotfire-test  at  Purdue  University  measured  ISP  and  C*  efficiency 
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2.2  TASK  II 

The  theoretical  performance  of  RTILs  can  be  superior  to  the  currently  used  fuels  for  both 
bipropellant  and  monopropellant  applications.  Implementation  by  the  Air  Force  of  such  RTILs 
would  be  in  accordance  with  the  green  propellant  initiative  and  would  meet  the  emerging 
demands  of  the  DoD.  While  RTILs  show  promise  as  enabling  energetic  propellants,  very  little  is 
understood  about  the  fundamental  reaction  chemistry  involved  in  the  energy  release  processes 
that  ensue  upon  decomposition/oxidation  of  these  fuels.  Achieving  reliable  ignition  and 
sustained  combustion  in  RTILs  are  important  challenges  for  the  Air  Force. 

In  this  effort,  we  have  carried  out  complementary  experimental  and  theoretical  studies  on  a 
select  few  model  RTILs  and  related  compounds  to  gain  fundamental  insight  into  the  reaction 
chemistry  of  these  substances.  Below,  we  provide  a  summary  of  our  findings.  Detailed 
descriptions  of  the  work  performed  and  the  conclusions  reached  may  be  found  in  articles 
published  during  the  tenure  of  this  research. 


TASK  II  Research  goals  and  objectives  FY14-FY16. 

In  order  to  model  the  ignition  behavior  of  energy  dense  materials,  such  as  room-temperature 
ionic  liquids  (RTILs),  an  important  goal  for  this  research  effort  was  to  improve  the  understanding 
of  the  reactive  intermediate  chemistry  involved.  The  following  objectives  were  pursued  towards 
this  end: 

1.  Elucidate  the  nature  of  the  reaction  mechanism  involved  in  the  thermolysis  of  model  RTILs. 

2.  Elucidate  the  nature  of  the  reaction  mechanism  involved  in  the  oxidation  of  DCA-based  RTILs. 

3.  Investigate  reaction  kinetics  and  flammability  envelopes  of  hypergolic  mixtures. 


Thermal  decomposition  of  CN  ionic  liquids.10 

The  thermal  decomposition  of  alkylimidazolium  ionic  liquids  with  cyano-functionalized  anions 
has  been  investigated  by  multiple,  complementary  experimental  techniques  and  density 
functional  theory  (DFT)  conductor-like  polarized  continuum  model  (CPCM)-generalized  ionic 
liquid  (GIL)  calculations.  Due  to  the  unusually  high  heats  of  vaporization  of  RTILs,  volatilization 
of  RTILs  through  thermal  decomposition  and  vaporization  of  the  decomposition  products  can  be 
significant.  Upon  heating  of  cyano-functionalized  anionic  RTILs  in  vacuum,  their  gaseous 
products  were  detected  experimentally  via  tunable  vacuum  ultraviolet  photoionization  mass 
spectrometry  performed  at  the  Chemical  Dynamics  Beamline  9.0.2  at  the  Advanced  Light  Source. 
Experimental  evidence  for  di-  and  tri-alkylimidazolium  cations  and  cyano-functionalized  anionic 
RTILs  confirms  thermal  decomposition  occurs  primarily  through  two  pathways:  (1)  deprotonation 
of  the  cation  by  the  anion  and  (2)  dealkylation  of  the  imidazolium  cation  by  the  anion  (Figure  20). 
Secondary  reactions  include  possible  cyclization  of  the  cation  and  C2-substitution  on  the 
imidazolium,  and  their  proposed  reaction  mechanisms  are  discussed  in  detail  in  Ref  10.  Two 
common  thermal  decomposition  mechanisms  have  been  confirmed,  carbene  formation  and  alkyl 
abstraction,  which  are  related  to  anion  basicity  and  nucleophilicity,  respectively.  Here  we 
propose  a  third  mechanism  that  likely  proceeds  through  the  carbene,  which  allows  for  addition 
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of  S  or  -NCN  groups  to  the  C2  position  on  the  imidazole  ring  and  could  be  useful  for  synthesis  of 
substituted  C2-imidazoles.  Additional  evidence  supporting  these  mechanisms  was  obtained 
using  thermal  gravimetric  analysis/mass  spectrometry,  gas  chromatography/mass  spectrometry 
and  temperature-jump  infrared  spectroscopy. 
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Figure  20.  Thermal  decomposition  mechanisms  of  EMI M+SCN',  including  -CH3  and  -CH2CH3 
abstractions  and  S  substitution  at  EMIM  C2  evidenced  by  vacuum  ultraviolet-time  of  flight  mass 
spectrometry. 

In  order  to  predict  the  overall  thermal  stability  in  these  ionic  liquids,  the  ability  to  accurately 
calculate  both  the  basicity  of  the  anions  and  their  nucleophilicity  in  the  ionic  liquid  is  critical.  Both 
gas-phase  and  condensed-phase  (CPCM-GIL)  density  functional  theory  calculations  support  the 
decomposition  mechanisms  and,  the  CPCM-GIL  model  could  provide  an  accurate  means  to 
determine  thermal  stabilities  for  ionic  liquids  in  general.  M06  density  functional  calculations  in 
the  gas  phase  and  the  SMD  (continuum  solvation  model)-GIL  variant  of  the  CPCM  for  the 
condensed  phase  have  demonstrated  the  ability  to  predict  trends  in  anion  basicity  and  anion 
nucleophilicity  in  pure  ionic  liquids,  which  should  prove  useful  in  predicting  thermal  stability 
trends  in  dialkylimidazolium  ionic  liquids  and  could  be  used  as  a  higher  accuracy  method  than 
the  gas-phase  DFT  approach  for  predicting  thermal  stabilities  of  ionic  liquids  in  general.  One 
important  finding  from  the  comparison  of  the  gas-phase  basicities  relative  to  the  GIL  condensed- 
phase  basicities  is  that  DCA'  is  more  basic  than  NO3'  in  the  condensed  phase  ionic  liquid,  which 
indicates  that  the  proton  transfer  from  HNO3  to  DCA~  is  likely  the  first  step  in  the  hypergolic 
ignition  mechanism  (Table  2). 
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Table  2.  Calculated  free  energies  of  acidity,  AGaCid,  in  the  gas  phase  and  in  the  condensed  phase 
by  SMD-GIL  and  SMD  (water)  at  the  M06/6-31+G(d,p)  level  of  theory.  Note  that  while  HDCA 
(HNCNCN)  is  more  acidic  than  HNO3  in  the  gas  phase,  in  the  SMD-GIL  condensed-phase,  HNO3  is 
more  acidic  than  HDCA,  possibly  facilitating  proton  transfer  from  HNO3  to  DCA':  HNO3  +  DCA-> 
NO3-  +  HDCA. 


AG  acid  (g) 

AG  acid  (1) 

AG  acid  (1) 

acid 

(kj/mol) 

SMD-GIL  (kj/mol) 

SMD-H2O  (kj/mol) 

HNCS 

1332.5 

HNCS 

607.0 

HNCS 

549.5 

HNO3 

1306.7 

HNCNCN 

578.0 

HNCNCN 

521.8 

HSCN 

1274.2 

HNO3 

566.6 

Hvdca  04 

515.2 

HNCNCN 

1272.6 

Hvdca  01 

565.2 

Hvdca  01 

514.6 

NCNHCN 

1234.6 

Hvdca  04 

564.9 

NCNHCN 

495.9 

Hvdca  01 

1212.1 

NCNHCN 

550.3 

HNO3 

494.9 

Hvdca  04 

1203.3 

HSCN 

548.8 

HSCN 

492.6 

HTCM  (central) 

1197.1 

HTCM  (central) 

537.9 

HTCM  (central) 

490.5 

HTCM  (terminal) 

1188.7 

HTCM  (terminal) 

533.9 

HTCM  (terminal) 

483.9 

Hvdca  02 

1152.7 

Hvdca  03 

518.1 

Hvdca  03 

471.4 

Hvdca  03 

1146.1 

Hvdca  02 

517.6 

Hvdca  02 

468.6 

vdca  =  vinylogous  dicyanamide,  Hvdca  01  =  terminal  NH,  Hvdca  02  =  -CCNiH,  Hvdca  03  =  -CCN2H,  Hvdca 
04  =  central  NH 


Reaction  of  aerosolized  DCA-ILs.11 

The  unusually  high  heats  of  vaporization  of  RTILs  complicate  the  utilization  of  thermal 
evaporation  to  study  ionic  liquid  reactivity.  Although  effusion  of  RTILs  into  a  reaction  flow-tube 
or  mass  spectrometer  is  possible,  competition  between  vaporization  and  thermal  decomposition 
of  the  RTIL  can  greatly  increase  the  complexity  of  the  observed  reaction  products.  In  order  to 
investigate  the  reaction  kinetics  of  a  hypergolic  RTIL,  l-butyl-3-methylimidazolium  dicyanamide 
(BMIM+DCA“)  was  aerosolized  and  reacted  with  gaseous  nitric  acid,  and  the  products  were 
monitored  via  tunable  vacuum  ultraviolet  photoionization  time-of-flight  mass  spectrometry  at 
the  Chemical  Dynamics  Beamline  9.0.2  at  the  Advanced  Light  Source.  Reaction  product 
formation  at  m/z  42,  43,  44,  67,  85, 126,  and  higher  masses  was  observed  as  a  function  of  HNO3 
exposure.  The  identities  of  the  product  species  were  assigned  to  the  masses  on  the  basis  of  their 
ionization  energies. 

To  our  knowledge,  this  is  the  first  study  on  the  chemical  reactivity  of  ionic  liquid  aerosols. 
Measurement  of  the  product  profiles  as  a  function  of  HNO3  exposure  allows  for  insight  into  the 
kinetics  of  the  initial  reaction  steps  in  the  hypergolic  reaction  of  BMIM+DCA~  +  HNO3.  Kinetics 
analysis  indicates  that  the  initial  reaction  of  BMIM+DCA~  with  HNO3  is  much  faster  than  the 
subsequent  reaction  of  HDCA  with  HNO3.  The  observed  exposure  profile  of  the  m/z  67  (HDCA+) 
signal  suggests  that  the  excess  gaseous  HNO3  initiates  rapid  reactions  near  the  surface  of  the  RTIL 
aerosol  (Figure  21).  Nonreactive  molecular  dynamics  simulations  support  this  observation, 
suggesting  that  diffusion  within  the  particle  may  be  a  limiting  step.  The  mechanism  is  consistent 
with  previous  reports  that  nitric  acid  forms  a  protonated  dicyanamide  species  in  the  first  step  of 
the  reaction,12  and  this  work  is  the  first  to  directly  detect  the  formation  of  HDCA.  Additional  FTIR 
experimental  results  indicate  that  NC>2+  appears  not  to  be  an  important  species  in  the  hypergolic 
ignition  mechanism. 
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Figure  21.  Depiction  of  the  diffusion-limited  process  in  the  ignition  of  hypergolic  DCA-based  ionic 
liquids  with  HN03- 

Reactivity  of  B-  and  Al-nanoparticle  infused  ionic  liquids.13 

The  interaction  of  H-functionalized  boron  nanoparticles  with  alkenes  and  nitrogen-rich  ILs  was 
investigated  by  a  combination  of  X-ray  photoelectron  spectroscopy,  FTIR  spectroscopy,  dynamic 
light  scattering,  thermogravimetric  analysis,  and  helium  ion  microscopy.  Surface  B-H  bonds  are 
shown  to  react  with  terminal  alkenes  to  produce  alkyl-functionalized  boron  particles.  The 
interaction  of  nitrogen-rich  ILs  with  the  particles  appears,  instead,  to  be  dominated  by 
boron-nitrogen  bonding,  even  for  ILs  with  terminal  alkene  functionality. 


Figure  22.  Ignition  delay  time  of  boron  nanoparticle-infused  ionic  liquids  MAT+DCA'  and 
AMIM+DCA'  as  a  function  of  particle  loading  as  determined  by  the  evolution  of  CO  2  by  rapid- 
scan  Fourier-transform  infrared  spectroscopy. 

This  chemistry  provides  a  convenient  approach  to  producing  and  capping  boron  nanoparticles 
with  a  protective  organic  layer,  which  is  shown  to  protect  the  particles  from  oxidation  during  air 
exposure.  By  controlling  the  capping  group,  particles  with  high  dispersibility  in  nonpolar  or  polar 
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liquids  can  be  produced.  For  the  particles  capped  with  ILs,  the  effect  of  particle  loading  on 
hypergolic  ignition  delaytimes  of  the  ILs  is  reported  in  Figure  22. 

Milling  of  boron  in  H2,  followed  by  mix-milling  with  a  capping  agent,  is  shown  to  be  an  efficient 
way  to  produce  boron  nanoparticles  capped  with  either  alkyl  groups  or  ionic  liquids.  IR 
spectroscopy  suggests  that  the  alkyl  capping  process  involves  reaction  of  terminal  C=C  bonds  in 
alkene  capping  agents  with  B-H  bonds  on  the  surface,  resulting  in  B-C  bond  formation.  For  the 
ionic  liquids,  the  mode  of  binding  to  the  surface  is  more  difficult  to  determine  unambiguously. 
There  is  no  obvious  difference  in  the  mode  of  binding  of  ILs  with  and  without  terminal  C=C  bonds 
(AMIM+DCA'  and  MAT+DCA',  respectively),  and  the  strong  perturbation  of  the  CeN  stretching 
vibrations  of  the  DCA"  anion  in  the  capping  layer  suggests  that  boron-DCA  interaction  is 
important  in  the  bonding.  Both  IR  and  zeta  potential  measurements  suggest  that  for  MAT+DCA' 
-capped  particles,  the  mode  of  binding  is  somewhat  different  for  particles  with,  and  without 
hydrogen  termination. 

Both  alkyl-  and  IL-capped  particles  are  protected  against  air  oxidation,  even  if  the  samples  are 
repeatedly  ultrasonicated  in  solvents  to  remove  any  weakly  complexed  capping  agents.  Several 
mechanisms  may  contribute  to  the  oxidation  resistance.  The  organic  capping  layers  may 
physically  exclude  O2  and  water  (the  primary  atmospheric  oxidants)  to  some  extent,  however, 
given  that  the  samples  were  heated  overnight  to  350  K,  (octene  boiling  point  =  394  K),  and 
exposed  to  air  for  at  least  6  hr  prior  to  XPS  analysis  (Figure  23),  it  seems  unlikely  the  capping  layer 
could  simply  be  functioning  as  a  diffusion  barrier.  More  likely,  bonding  of  the  capping  layer  to 
the  surface  changes  the  surface  chemistry,  such  that  it  is  unreactive  to  O2  and  water  up  to  350  K. 


Figure  23.  X-ray  photoelectron  spectrum  (XPS,  left)  of  the  oxidized  and  protected  boron 
nanoparticles,  and  scanning  electron  microscope  image  (SEM,  right)  of  the  ball  milled  boron 
nanoparticles  that  have  been  protected. 
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MD  calculations  of  thermal  decomposition  of  DNB.14 

Certain  room-temperature  ionic  liquids  exhibit  hypergolic  activity  as  liquid  bipropellants. 
Understanding  the  chemical  pathways  and  reaction  mechanisms  associated  with  hypergolic 
ignition  is  important  for  designing  new  fuels.  It  has  been  proposed12  that  an  important  ignition 
step  for  the  hypergolic  ionic  liquid  bipropellant  system  of  dicyanamide/nitric  acid  is  the  activation 
and  dissociation  of  the  1,5-dinitrobiuret  anion  DNB'.  Forthe  work  reported  here,  a  quasiclassical 
direct  dynamics  simulation,  at  the  DFT/M05-2X  level  of  theory,  was  performed  to  model  H+  + 
DNB'  association  and  the  ensuing  unimolecular  decomposition  of  HDNB.  This  association  step  is 
324  kcal/mol  exothermic,  and  the  most  probable  collision  event  is  for  H+  to  directly  scatter  off  of 
DNB',  without  sufficient  energy  transfer  to  DNB'  for  H+  to  associate  and  form  a  highly 
vibrationally  excited  HDNB  molecule.  However,  approximately  1/3  of  the  trajectories  do  form 
HDNB,  which  decomposes  by  eight  different  reaction  paths  and  whose  unimolecular  dynamics  is 
highly  nonstatistical.  Some  of  these  paths  are  the  same  as  those  found  in  a  direct  dynamics 
simulations  study  of  the  high-temperature  thermal  decomposition  of  HDNB15  for  a  similar  total 
energy. 

Table  3.  Direct  dynamics  trajectory  simulation  results  for  H*  +  DNB*  collisions. 
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The  reaction  pathways  found  in  the  direct  dynamics  simulation  modeling  of  H+  +  DNB~ 
collisions  are  shown  in  Table  3.  Possible  formation  of  the  highly  vibrationally  excited  HDNB 
molecule  is  modeled  by  assuming  an  electronic  transition  from  the  H+  +  DNB'  to  H*  +  *DNB 
potential  energy  surface  and  simulating  the  collisions  on  this  latter  surface.  The  most  likely 
event,  that  is,  2/3  of  the  trajectories,  is  for  the  H-atom  to  directly  scatter  off  of  »DNB  without 
forming  a  "hot"  HDNB  molecule.  For  most  of  the  trajectories,  an  insufficient  amount  of  the  H* 
+  *DNB  relative  translational  energy  is  transferred  to  »DNB  vibrational  energy  to  form  a 
vibrationally  excited  HDNB  molecule.  Another  way  to  view  these  dynamics  is  that  there  is 
insufficient  IVR  (internal  vibrational  energy  redistribution)  from  the  incipient  H-DNB  bond  of 
the  H*  +  «DNB  collision  to  other  vibrational  modes,  and  the  H«  atom  directly  scatters  from 

•  DNB.  Approximately  1/4  of  the  trajectories  followed  path  2  in  Table  3,  forming  the  H2N2O2 
product,  which  has  an  isomerization  pathway  to  HNNOOH.  The  H2N2O2/  HNNOOH  product 
pathway  may  undergo  seven  different  secondary  dissociations.  This  reaction  occurs  by  the 
colliding  H-atom  first  interacting  with  the  central  N-atom  and  the  two  adjacent  O-atoms  of 

•  DNB,  which  is  then  followed  by  IVR  promoting  the  chemical  reaction.  This  was  found  to  also 
be  an  important  pathway  for  HDNB  dissociation  in  a  previous  direct  dynamics  simulation  by  Liu 
et  al.,10b  for  which  HDNB  was  thermally  excited.  In  summary,  a  comparison  of  the  current 
results  with  the  previous  study  by  Liu  et  al.10b  shows  that  the  initial  conditions  for  the  excited 
HDNB  moiety  are  crucial  for  its  ensuing  unimolecular  dynamics.  H+  +  DNB'  association  localizes 
the  energy  in  HDNB,  creating  a  "hot  spot",  and  IVR  is  not  sufficiently  rapid  to  result  in  the 
unimolecular  dynamics  found  by  Liu  et  al.  for  their  random,  thermal  initial  conditions.  For  this 
previous  simulation,  there  were  rapid  statistical  interconversions  between  all  of  the  HDNB 
conformers  before  unimolecular  dissociation  occurred.  In  future  work,  it  would  be  of  interest 
to  study  both  the  random  and  nonrandom  excitation  of  HDNB  in  a  condensed  phase  liquid 
environment. 

High-level  a b  initio  calculations  of  DNB.16 

Mononitrobiuret  (MNB)  and  1,5-dinitrobiuret  (DNB)  are  tetrazole-free,  nitrogen-rich,  energetic 
compounds.  For  the  first  time,  a  comprehensive  ab  initio  kinetics  study  on  the  thermal 
decomposition  mechanisms  of  MNB  and  DNB  has  been  carried  out.  In  particular,  the 
intramolecular  interactions  of  amine  H-atom  with  electronegative  nitro  O-atom  and  carbonyl  O- 
atom  have  been  analyzed  for  biuret,  MNB,  and  DNB  at  the  M06-2X/aug-cc-pVTZ  level  of  theory. 
The  results  show  that  the  MNB  and  DNB  molecules  are  stabilized  through  six-member-ring 
moieties  via  intramolecular  H-bonding  with  interatomic  distances  between  1.8  and  2.0  A,  due  to 
electrostatic  as  well  as  polarization  and  dispersion  interactions.  Furthermore,  it  was  found  that 
the  stable  molecules  in  the  solid  state  have  the  smallest  dipole  moment  amongst  all  the 
conformers  in  the  nitrobiuret  series  of  compounds,  thus  revealing  a  simple  way  for  evaluating 
reactivity  of  fuel  conformers. 

The  PESs  for  thermal  decomposition  of  MNB  and  DNB  were  characterized  at  the  M06-2X/aug- 
cc-pVTZ  and  RCCSD(T)/cc-pV°°Z//M06-2X/aug-cc-pVTZ  level  of  theory.  In  particular,  the  values 
of  the  energy  barriers  and  endothermicities  at  the  M06-2X/aug-cc-pVTZ  level  of  theory  show 
remarkable  agreement  with  the  values  obtained  from  the  RCCSD(T)/cc-pV°°Z//M06-2X/aug-cc- 
pVTZ  computations,  implying  that  the  former  level  of  theory  could  be  applicable  to  larger 
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analogous  systems  (Figure  24).  It  was  found  that  the  thermal  decomposition  of  MNB  is  initiated 
by  the  elimination  of  HNCO  and  HNN(0)0H  intermediates.  Intramolecular  transfer  of  a  H-atom, 
respectively,  from  the  terminal  NH2  group  to  the  adjacent  carbonyl  O-atom  via  a  six-member- 
ring  transition  state  eliminates  HNCO  with  an  energy  barrier  of  35  kcal/mol  and  from  the  central 
NH  group  to  the  adjacent  nitro  O-atom  eliminates  HNN(0)0H  with  an  energy  barrier  of  34 
kcal/mol.  Elimination  of  HNN(0)0H  is  also  the  primary  process  involved  in  the  thermal 
decomposition  of  DNB,  which  processes  C2v  symmetry.  The  energy  barrier  for  HNN(0)0H 
elimination  in  DNB  is  6.60  kcal/mol  lower  than  that  in  MNB  due  to  an  extra  hydrogen  bond  in  the 
transition  state  for  the  former,  which  results  in  DNB  being  less  stable  than  MNB.  Furthermore, 
the  HNN(0)0H  intermediate  subsequently  decomposes  via  multiple  wells  and  multiple  channels. 
The  RRKM/multi-well  master  equation  simulations  revealed  that  decomposition  to 
thermodynamically  stable  N2O  +  H2O  products  via  isomerization  is  the  primary  channel  during 
HNN(0)0H  decomposition.  The  rate  coefficients  for  the  primary  decomposition  channels  for 
MNB  and  DNB  were  quantified  as  functions  of  temperature  and  pressure. 

The  rate  coefficient  data  can  be  used  to  interpret  the  experimental  results  of  Klapotke  et  al.17 
regarding  the  thermal  stability  of  MNB  and  DNB,  and  their  decomposition  products  and  ignition. 
Such  information  is  essential  in  the  design  and  manipulation  of  molecular  systems  for  the 
development  of  new  energetic  materials  for  advanced  propulsion  applications. 


Figure  24.  The  potential  energy  surface  for  the  thermal  decomposition  of  DNB  calculated  at  the 
M06-2X/aug-cc-pVTZ  level  of  theory. 


N2H3  +  NO2  reaction  kinetics.18 

The  N2H3  +  NO2  reaction  plays  a  key  role  during  the  early  stages  of  hypergolic  ignition  between 
N2H4  and  N2O4.  Here  for  the  first  time,  the  reaction  kinetics  of  N2H3  in  excess  NO2  was  studied  in 
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2.0  Torr  of  N2  and  in  the  temperature  range  298-348  K  in  a  pulsed  photolysis  flow-tube  reactor 
coupled  to  a  mass  spectrometer.  The  temporal  profile  of  the  product,  HONO,  was  determined 
by  direct  detection  of  the  m/z  +47  ion  signal.  For  each  chosen  [NO2],  the  observed  [HONO]  trace 
(Figure  25)  could  be  fitted  to  a  bi-exponential  kinetics  expression,  which  yielded  a  value  for  the 
pseudo-first-order  rate  coefficient,  k',  for  the  reaction  of  N2H3  with  NO2.  The  slope  of  a  plot  of  k' 
versus  [NO2]  yielded  a  value  for  the  total  bimolecular  rate  coefficient,  which  could  be  fitted  to  an 
Arrhenius  expression;  (2.36  ±  0.47)  x  10'12exp((520  ±  350)/T)  cm3  molecule  1  s'1.  The  errors  are 
la  and  include  estimated  uncertainties  in  the  NO2  concentration. 


Time  (sec) 

Figure  25.  Typical  [HONO]  temporal  profile  observed  in  the  flow-tube  reactor  at  298  K. 

The  potential  energy  surface  of  N2H3  +  NChwas  investigated  by  several  advanced  ab  initio 
theories,  including  coupled-cluster  and  multi-reference  second-order  perturbation  methods,  and 
the  results  reveal  a  new  reaction  mechanism.  The  reaction  is  exothermic  by  up  to  42  kcal/mol, 
and  proceeds  via  a  complex  addition-isomerization-dissociation  mechanism  with  the  transition 
state  energies  of  all  the  reaction  channels  either  below  or  nearly  equal  to  that  of  the  entrance 
asymptote,  implying  all  reaction  channels  are  important  for  products  formation. 

It  was  found  that  the  direct  addition  of  NO2  to  one  side  of  N2H3,  via  a  6-member-ring  transition 
state,  forms  the  adduct,  N2H3NO2,  and  that  addition  to  the  other  side  of  N2H3,  via  a  5-member¬ 
ring  transition  state,  forms  a  complex  that  undergoes  two  facile  isomerization  reactions  to  form 
the  adduct,  N2H3ONO.  The  subsequent  isomerization  and  decomposition  of  the  addition  adducts 
lead  to  various  products.  The  temperature  and  pressure  dependent  rate  coefficients  for  all  of 
the  isomerization  and  dissociation  channels  were  computed  via  RRKM/multi-well  master 
equation  simulations  (Figure  26).  The  results  reveal  that  the  dominant  channels  for  the  N2H3  + 
NO2  reaction  system  lead  to  the  formation  of  the  products,  trans-HONO  +  trans-NH=NH,  and  NO 
+  NH2NHO,  through  the  addition  adducts  N2H3NO2  and  N2H3ONO,  and  these  addition  adducts 
show  very  strong  negative  temperature  dependences  in  the  temperature  range  from  300  to 
3000K. 
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Figure  26.  High  pressure  limit  rate  coefficients  for  the  dissociation  and  isomerization  channels 
of  N2H3  +  NO2  —>  N2H3NO2  —>  Products  (left)  and  of  N2H3  +  NO2  —>  N2H3ONO  —>  Products  (right). 

Flammability  limits.19 

A  numerical  method  is  demonstrated  in  which  a  simple  flame  temperature  criterion  of  2700  K 
is  used  to  map  out  flammability  diagrams  as  a  function  of  total  mixture  pressure  and  equivalence 
ratio  in  the  hypergolic  system,  monomethylhydrazine/nitrogen  tetroxide/helium 
(MMH/NTO/He).  The  computed  results  are  in  good  agreement  with  experimentally  determined 
ignition  diagrams  for  MMH/NTO/He  (Figure  27).  The  method  is  used  to  predict  the  lower  and 
upper  flammability  limits  of  other  hypergolic  mixtures  at  298  Kand  1  atm.  This  numerical  method 
allows  for  the  rapid  assessment  of  the  flammability  hazards  of  uncharacterized  fuel/oxidant 
mixtures  that  may  be  encountered  in  the  work  place  as  well  as  their  auto-ignitability 
(hypergolicity)  in  combustion  devices. 


Figure  27.  MMH/NTO/He  (with  60  mol%  of  He)  flammability  diagram  as  a  function  of  total  mixture 
pressure  and  equivalence  ratio  obtained  by  assuming  rapid  mixing.  EH :  no  ignition;  ^ :  ignition. 
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Future  work  will  extend  this  method  to  mixtures  of  ionic  liquid  fuels  with  common  oxidizers 
such  as  HNO3,  NTO,  and  H2O2.  This  method  can  be  utilized  to  systematically  screen  for  suitable 
ionic  liquid/oxidizer  candidates  for  potential  hypergolic  propellant  applications  provided  that  the 
chemical  kinetics  model  for  ignition  is  sufficiently  developed  to  determine  temperature  profiles 
to  assess  that  "reasonable"  i.e.,  not  too  long,  ignition  delays  are  expectable. 


Conclusion 

In  this  work  we  have  used  a  variety  of  complementary  experimental  techniques  and  theoretical 
approaches  to  elucidate  the  reaction  mechanisms  involved  in  the  decomposition  of  energetic 
RTILs  as  a  result  of  thermolysis,  catalysis  and  oxidation.  As  our  knowledge  improves  in  these 
areas,  it  should  be  possible  to  build  predictive  numerical  models  for  the  accurate  assessment  of 
the  performance  and  the  state-of-health  in  RTIL  monopropellant  and  bipropellant  thrusters. 
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USAF-AFRL 

Brand 

NASA 

Providing  advanced 
monopropellant  based  on 
AFOSR  material  for  spacecraft 
demonstration 

Spacecraft 

Monopropulsion 

L,l 

O 

Pd 

IONIC  LIQUID- 
BASED 

PROPELLANTS 

6/2012- 

12/2014 

BERMAN 
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1 
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1 
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Transitioned  From:  AFRL  =  L;  Industry  =  I;  Academia  =  A 

Transitioned  To:  Industry  =  I;  Air  Force  6.2  or  6.3  =  AF;  Other  AF,  DoD,  Government,  etc.  =  O  (please  specify) 

Application:  Product  (New  or  Improved)  =Pd;  Process  (New  or  Improved)  =  Pc;  Other  =  O  (please  specify) 


Distribution  A:  Approved  for  Public  Release;  Distribution  unlimited 
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